Carrier mediated interlayer exchange, ground state phase diagrams and transition 

temperatures of magnetic thin films 
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We investigate the influence of the carrier density and other parameters on the interlayer exchange 
in magnetic thin film systems. The system consists of ferromagnetic and non-magnetic layers where 
the carriers are allowed to move from layer to layer. For the ferromagnetic layers we use the Kondo- 
lattice model to describe interactions between itinerant electrons and local moments. The electrons' 
properties are calculated by a Green's function's equation of motion approach while the magnetiza- 
tion of the local moments is determined by a minimization of the free energy. As results we present 
magnetic phase diagrams and the interlayer exchange over a broad parameter range. Additionally 
we can calculate the transition temperatures for diflFerent alignments of the ferromagnetic layers' 
magnetizations . 



I. INTRODUCTION 



Magnetic thin films received much attention in 
experimentalii^ and theoretical^i^ research, especially af- 
ter the discovery of a giant magneto resistance (GMR^i^). 
One result of this discovery was a massive increase of 
storage density in hard drives. The phenomenon of the 
GMR is based on different possible alignments of mag- 
netizations in the single layers. If they are parallel to 
each other the resistance for a current perpendicular to 
the layer is considerably smaller than in the antiparal- 
lel alignment. Thus it is very interesting to know under 
which conditions a film material favors to align its layer 
magnetizations parallel or antiparallel. 
The way the magnetizations are aligned is defined by the 
so called interlayer exchange, which can either be ferro- or 
antiferromagnetic. Even though the interlayer exchange 
in most materials is assumed to be carrier-mediated^'^ it 
is usually described by a Heisenberg model via an RKKY 
interaction^"^^. This model successfully describes the 
damped oscillations of the interlayer exchange accord- 
ing to the number of layers^, which essentially influence 
the type of alignment of the magnetizations. 
On the other hand the standard RKKY interaction is 
derived from the free electron gas, which underestimates 
the interactions in strongly correlated materials. Even 
modifications of the RKKY—, which incorporate these 
interactions much more, are based on the Heisenberg 
model. However, many materials are not best described 
by a Heisenberg model, but rather, e.g., by a Kondo- 
lattice/double exchange model. Famous examples are 
the magnetic 4/-systems (Gd, Tb,. . . ) or the mangan- 
ites (Lai_a:Caa:Mn03, Lai-^rSra^MnOa,. . . ). Contrary to 
the Heisenberg model the Kondo-lattice model (KLM) 
directly describes itinerant carriers aside the local mo- 
ments. Because the itinerant carriers play an important 
role for many materials we do not want to use the RKKY 
method to map the KLM onto a Heisenberg model, since 
it is unsure how much the original properties of the KLM, 
and with it the influence of the carriers, are changed 
within this procedure. 



The KLM consists of two subsystems - itinerant electrons 
and local moments. It is difficult to treat both subsys- 
tems in one theory on the same level which is one reason 
for the mapping onto a Heisenberg model in the RKKY 
method to determine the local moments' properties. In 
our method, similarly to the RKKY method, we first cal- 
culate the electronic properties while we leave the local 
moments' values (e.g. magnetization {S^)) as parame- 
ters. But contrary to the RKKY we then calculate the 
free energy of the system as a function of the magnetiza- 
tion of the local moments^^ and find its minimum. Thus 
we need no mapping on the Heisenberg model. 
Generally the KLM describes an on-site interaction of lo- 
cal moments and an itinerant electron's spin cr^. De- 
pending on the coupling J they prefer to align either par- 
allel (J > 0) or antiparallel (J < 0). Due to the hopping 
of the electrons the whole system can become magneti- 
cally ordered. The Hamiltonian of the KLM reads 



(1) 



where Tij is the hopping integral and the c^^^ are the 
creators (annihilators) of electrons with spin a — ±1 at 
the lattice site Ri. It is convenient to write the scalar 
product Si ■ (Ti explicitly by the use of the Pauli matri- 
ces' vector or and cri = \ cf^^aa'Cicr' to get the 
second line of ([T]). There we additionally introduced 
Za — — 6a, -1, riic — cf^Cia and the spin ladder 

operator' S'f = Sf + iz^Sf. 

One of the main purposes of this work is to calculate the 
interlayer exchange in dependence on important model 
parameters. We will only treat "sandwiched" layers in 
this paper. In particular the system consists of Nl layers 
where the top and the bottom layers are ferromagnetic 
(FM), i.e. described by the KLM (cf. Fig. [T]). These 
two FM layers are separated from each other by Nl — 2 
non-magnetic (NM) layers (spacers). The FM layers can 
"communicate" with each other by means of the inter- 
layer exchange. 
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(Nl-2) NM 
layers 




Figure 1: Schematic picture of tiie layer system investigated 
in this work. The system consists of Nl layers with two fer- 
romagnetic (FM) layers at the top and the bottom. They are 
separated by Nl — 2 non-magnetic (NM) layers. 



Due to the interlayer exchange the magnetizations of 
both FM layers can align in different ways. We restrict 
the types of alignment to be either parallel or antiparal- 
lel. Other phases like canted ones are neglected in our 
work. Since the layer system is symmetric the magneti- 
zations should be equal in magnitude, which reduces the 
amount of numerical calculations significantly. 



A. Electronic subsystem 

We want to find an approximative solution for the elec- 
tron Green's function 



G.,.(i?)=((c..;c+)) 



or its Fourier-transformed equivalent 

Gkcr 



{E) 



(4) 



(5) 



To do this we calculate several equations of motion which 
reads for a general Green's function with the operators 
A and B 



E{{A-B)) = {[A,B]^) + {{[A,'HUB)) 



(6) 



where (...) is the thermodynamic average. The com- 
mutator [A,T-C\_ on the right-hand-side usually leads to 
higher Green's functions which include more operators 
and prevent a direct solution of the problem. Thus we 
have to find a physically reasonable decoupling scheme to 
trace back these higher functions to those of lower order. 



1. Film systems 



II. THEORY 



To compare different phases a system can exhibit the 
key quantity is the internal energy. In the Kondo-lattice 
model it is simply given by 



U 
N 



^ 1 di? Ef^{E)p„{E) 



(2) 



where f-{E) is the Fermi function and Pa{E) the quasi 
particle density of states (QDOS). The QDOS can be 
calculated from the electronic Green's function Gkcr(-E) 
via the condition 



r{E) 



(3) 



Thus we need to find an expression for the Green's func- 
tion, which will in general depend on electronic param- 
eters (e.g. band occupation) and on the spins' proper- 
ties (e.g. magnetization). First we want to calculate the 
Green's function only in the electronic subsystem and 
treat spin expectation values as parameters. In a second 
step we will show how the optimal spin expectation val- 
ues can be determined by a minimization of the system's 
free energy. 



First we want to discuss the general procedure in a film 
systemic. The Hamiltonian can be written as 



a/3 + 



lA 



(7) 



with a general interaction term H^^ . The lattice indices 
i II , j|| refer to lattice sites Ri^ within a film plane while 
a, P name a distinct layer. Different layer positions are 
described by a primitive translation r" perpendicular to 
the film planes (R°J| = Ri|| -I- r"). A movement of the 

electrons is incorporated by a hopping term ^i^^n where 
the electron can move within the film from Rf to Rf 



or between two layers a, B (Rf — )■ Rf = Rf + r" 

J ^ I A" V l|| ^ I|| i|| 

To simplify the calculations we restrict the movement to 
next-neighbor hopping in a simple cubic lattice 



»iiJii ~ 8 «iiJii±i 



W 



2D 



tutu 



(8) 



with S, 



'iiiJii ^^^^ band width of a single film 

W^^ and the resulting dispersion 



1 



rpal3 ik||(Ri|| 



(9) 



The Fourier transformation is restricted to the layer with 
iV|| sites and the 2-dimensional wave vector k|| . We can 
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set up the equation of motion for the electron Green's 
function 



T, 



7/3 



G 



7/3 



(10) 
(11) 



'llJll 



lAl 



In the last step we introduced a self-energy Maa{E) 
which leads to a formal solution of the equation. Gen- 
erally the self-energy can be non-local, ^i'^^'j^^aiE), but 
only local ones will occur in this work. With a Fourier- 
transformation as in we find the Green's function by 
an inversion of an (N^ x iVi)-matrix 



(12) 



9^a^{E) 



with 



9^U{E) = E - - MUE) , 



...J 



(13) 
(14) 



When we find an expression for the self-energy the prob- 
lem is solved. 



2. Kondo lattice system 

To find a solution for the self-energy, we explicitly use 
the Kondo-lattice model adapted to film geometry now: 



(15) 



a/3 
J 

~ 2 



Up to now no general solution for the KLM is known 
and that is why we use a moment conserving decou- 
pling approach (MCDA^^) to get the self-energy Macr{E). 
This non-perturbative approach is based on equations of 
motions up to the second order and several decoupling 
schemes. It has been used successfully various times to 
describe different aspects of the KLMi^ii^iilii^. We want 
to sketch the general procedure and refer to Ref for 
details. Starting point is the equation of motion of the 



electron Green's function 

17 rp-t(3 



E 

ill ,7 



2 



cap _^f^ j-aa^ (' P^l 4_ p""'' f p\ 



(16) 

(17) 



with 



f: 



JE) = {{S-lo 



iiia^'"||7,-o-j 



(18) 
(19) 



Two new higher Green's functions appear; the Ising func- 
tion /,,"];f||^|^^(i;) and the spin-flip function i^^^^^^- ^^(i;). 

Due to the general expression of G"^j^^^{E) via a self- 
energy as in pT|) they have to be connected to Maa{E) 
as indicated in ([TT]). Thus when we find appropriate ex- 
pressions for both functions we get the self-energy, too. 
To do this we set up the equations of motion for both 
functions: 



E 



{es^:^ 



Tl 



'llJll 



[E] 



(20) 



Nat/^ I // 

'-'j||ay"^m||j|l ^ \\ 



'^i«a'^rn«'ya, 



lA 



(21) 



lA 



The commutators on the right-hand-side prevent a direct 
solution. We now differentiate between diagonal — 
m\\ and a = 7) and non-diagonal terms ^ m\\ or 
a ^ 7). For the non-diagonal ones we approximate 



{{Slab 



my 7cr ; 



H 



« M^„{E)l^ 



(22) 



ilimiiii 



AE) 



(23) 



Af,,_.(i?)F";;^,_,„,(i?) 



and 



St.. a ■> El 



lA 



Cm ||7(t; Cjii^ct/ 



lA 



« 
= 



In ((22)) . p3|) we used the formal equivalence 

[cm||7,±(7; ]_ ^ -^'-^7,±(7 (-^)c?n II 7±o 



(24) 
(25) 

(26) 



due to the Dyson equation, while the Green's functions 
in (1^^ . (psj) can be connected to magnon energies, which 
are much lower than the electronic energies^®. This holds 
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for all spin quantum numbers S. 

Due to higher correlations the same procedure cannot 
be done for the diagonal terms. But in this case we can 
calculate the commutators in ([20)) and (I?!]) explicitly. As 
usual new higher Green's functions appear 



i^}.(i?) = ((5f5fc..;c+)) 



(27) 
(28) 
(29) 
(30) 



where we suppressed the layer indices for readability 
{i = i||a, j = The first two functions can be writ- 

ten in the limiting cases of low spin S* = ^ or ferromag- 
netic saturation {S^} = S as linear combinations of the 
electron Green's function, the Ising function and/or the 
spin-flip function. Thus we assume that a linear com- 
bination is also possible aside these limiting cases and 
choose 



(1) 



(31) 
(32) 



The coefficients can be explicitly determined by the con- 
dition that the moments of these functions have to be 



conserved. Contrary to Refji^ we set 

i^]^,(£;)=((5>,.Q^;c+)) 



(33) 
(34) 



which effectively means a neglection of double occupa- 
tion effects, for the sake of simplicity. Recent results^ 
show that double occupation effects play a minor role 
at high spins. On the other hand even with inclusion 
of double occupancies, which is essential for lower spins, 
results (phase diagrams, transition temperatures) hardly 
depend on the spin quantum number S but on the cou- 
pling strength JS* for J > 0. Thus we will only use S" 3> ^ 
and J > later on, which allows the neglection of dou- 
ble occupancies while the dependence on JS holds for all 
spins. 

Due to the simplifications all terms are traced back onto 
four unknown functions: Gf^- „, F°"f^- „, iT'^^n „ and 

^»||J||Cr' HH\->\f' /il'll^ll'^ 

Maa ■ The equations of motion of the spin- flip and Ising 
function reduce to analytical expression which have the 
solutions 

X = I^F. (35) 

Both are proportional to the electron Green's function 
G°||'^||^(i?). The prefactors Xarj{E) only contain lo- 
cal terms, like expectation values, the local self-energy 
Maa{E) and the local Green's function 



With ([TT)) we can write the self-energy 

Mc^^iE) = {zJ^^{E) -t- F^AE)) (36) 

which is then indeed local, too. In fact it was not neces- 
sary to assume a local self-energy in (jl7p but it simplifies 
the derivation of Maa{E). The reason for the locality of 
the self-energy is only due to the neglection of magnon 
energies^. 

Especially due to the dependence of Maa {E) on the local 
electronic Green's function 

MaAE) = Maa{Ga,±AE)] ■ (37) 

the self-energy is through ([T2|) connected to the Green's 
functions of all other layers. Since we get all Green's 
functions Gg"^(i<^) from ([T^ the whole system of equa- 
tions is closed. AU implicitly included, e.g. in Oal 
in (|3ip . electronic expectation values {{naa), {Sanaa), 
. . . ) are obtainable by the spectral theorem from the 
involved Green's functions. The spins' expectation val- 
ues {{{S^Y){lp), ((S'^)^)(v9), ...) can be derived from 
the magnon density which wc can get from a given 
magnetizatio n^^'^" (5^) = M{ip) ^ fiM), which we see 
as a parameter. 

With the calculation of the Green's function we also get 
the important quantity of the internal energy 

^ ^ E / dE EUE)lmGZ;aiE) (38) 



N 



which can be used to compare the energies of different 
phases. 

In the special systems we treat in this work (cf. Fig. [1]) 
only the layers at the top and the bottom have a non- 
vanishing (direct) self-energy {Maa{E) = 0, Va 7^ 1, iV^), 
but also the NM layers are influenced by interactions 
within the FM layers due to a hybridization induced by 
the hopping in ([12]). Additionally the FM and NM 



layers can differ in their centers of gravity r""| of the 
free bands. Throughout the paper we use 



while 



rpaa °'i^hJ^L NM 



(39) 



(40) 



may be unequal zero. 



B. Local moments 

In the previous section we derived an expression for 
the electronic Green's function and thereby also for the 
internal energy. One drawback of the method is that we 
have to see the magnetization of the local moments as a 
parameter. Thus we cannot determine M = (S^) from 
the equations itself. This can be solved if we can get the 
free energy F{M,T) for many values of M and decide 
which magnetization leads to the lowest energy. 
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1. Free energy 

The free energy is given by 

F{T, M) = U{M, T) - TS{M, T) (41) 



with the entropy S{M,T). We can use Sm{T) 



— g^M- and get 



Fm{T)^Um{T)+T- 



dT 



which can be transformed via the product rule and 
^ _ dFM{T) 



^va{-{FM{T)-FM{T = Q)) 



= -SMiT = 0) 



(42) 



(43) 



to 



Fm{T) 



--Um{0) 



dT 



TSMiO)- 
,Um{T')^Um{0) 



(44) 



(TO 



/\2 



=Im{T) 

Since we can get the internal energy from the method 
of Sec. Ill Al the free energy can be calculated when we 
find an expression for the entropy at T = 0. When we 
assume this can be done we can calculate the free energy 
for many values of M for a distinct temperature T. Thus 
we get an ensemble of energies {Fm{T)} which leads to 
a magnetization dependent free energy 



{Fm{T)} ^ Ft{M) . 
When we find its minimum according to M 



dFriM') 



dM' 



= 



(45) 



(46) 



M 



we get the optimal magnetization of the system. 
Normally one would use the system's total magnetization 
M*°* = M + a^{M) of local moment and electron mag- 
netization cr^ for the minimization of the free energy, but 
since the total magnetization in our case is a function of 
M it is sufficient to use M as the order parameter. 



2. Entropy 

The aim of this section is to find an approximate ex- 
pression for the entropy at T = 0. Generally the entropy 
can be written as 



SoiM) = fcslnr 



M 



(47) 



where Tm is the number of states which contribute to 
a magnetization M and ks is the Boltzmann constant. 
Note that the expression (|47l) is not the ground state 



entropy! Due to the third law of thermodynamics the 
ground state entropy would be a constant according 
to the system parameters. But the system is not in 
the ground state for most magnetizations, but only for 
one optimal M°p*. Thus the ground state entropy is 
S'o(Af°P*) while the complete magnetization dependent 
entropy 5*0 ( Af) is not visible in reality but plays the role 
of an auxiliary quantity to find the optimal value M°p'. 
The whole systems consists of itinerant electrons and 
local moments which both affect the entropy. First 
we make the ansatz that the total number of states 
can be written as a product of the subsystems' states 
Tm — ^^m^'mi which means that the number of states 
of one subsystem is not affected by a special state of the 
other one. This ansatz can be motivated by the dynamics 
of the subsystems, which differ in several orders of mag- 
nitude between the fast electrons and the slow moments. 
Thus the total entropy is given by a sum of the single 
ones: 

So{M)^S'°'{M) + S'o\M) . (48) 

For the electron entropy we use the expression for a Fermi 
gas 

SfiM)^~kB J2 [(l-("k||«.))ln(l-(nk||,,))+ 

k|| ,a,cr 

+ ("-kiiacr) ln(nk||Qcr)] (49) 

which seems plausible since the electrons see more or 
less a static potential of the moments. The mean val- 
ues (n kjicto-) are, however, calculated by the full theory 
of Sec. ll^at a fixed M! 

Each local moment can have a spin projection mf = 
—S,—S + while the complete system can be 

described by a configuration {mf}. We now assume 
that all configurations with the same magnetization 
{mf\N~^ J2i = M} have the same energy while con- 
figurations with different magnetizations have different 
energies. This can be for example, but not exclusively, 
the case if the electrons act as an effective field on 
the local moments (ideal paramagnet). The size or sign 
of the field plays no role as long as B'^^ 7^ 0. It means 
that a special magnetization is equivalent to a distinct 
energy and we can count the number of states as in the 
microcanonical ensemble 

E<Ei<E+AE 

^ T^°%M) = J2 1 ■ (51) 

Thus the number of states is equal to the number of con- 
figurations which leads to a magnetization M. These as- 
sumptions can only be made because we want to find the 
entropy at T = 0. At larger temperatures more energy 
levels contribute to a thermodynamical mean value like 
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-2-1012 

E[eV] 

Figure 2: Spin up/down quasi particle density of states of 
independent layers for an FM {solid line, M = S) and an NM 
layer {dashed line). Parameters: S = 7.5, JS — 2eV, n = 0.5, 



lev, 







the magnetization due to energy/temperature dependent 
factors like - e'''^^)"'^'. 

We find that, for example, the entropy for spins with 
5 = i is given as 

So{M) 



;i-A^)In(i-M)^ 

+ (i + M)ln(i + M) 



(52) 



Note that in the film systems treated in this work not 
every lattice site has to be occupied by a local moment 
which reduces the local moment entropy. The complete 
lattice consists of = N^N^^ lattice site (A^y is the num- 
ber of sites of one layer). Due to the independence of 
local moments the total entropy can be expressed by a 
sum of layer entropies Sq'^°'^{M) 



clot 



(M) 



N 



(53) 



where the entropies of the NM layers are zero. 
The approximations used in this section for the entropy 
are reasonable for the Kondo lattice model, but cannot 
be assigned to other systems, e.g. the Heisenberg model, 
without further considerations of their special properties. 



III. RESULTS 

We want to discuss the existence of magnetic configu- 
rations. At vanishing temperature T = we only have to 
compare the internal energies which can be directly cal- 
culated from the electronic Green's function. At higher 
temperatures T > there exists an interplay between a 



reduction of internal energy and an increase of entropy. 
Therefore we have to use the free energy as the decisive 
thermodynamic potential in this case. 



A. Phase diagrams at T=0 

Before we treat coupled layer systems it is instructive 
to mention some properties of independent layers, since 
important characteristics are valid for the coupled ones, 
too. First of all it is a basic question whether or not 
a layer described by a Kondo-lattice model is actually 
ferromagnetic for a given parameter constellation. We 
want to concentrate on two parameters: the coupling J 
(or better JS) and the band occupation n. In the Kondo 
lattice model (with neglection of double occupancies) the 
electron band splits into two subbands which are for high 
spins at energies E w ^^JS (cf. Fig. [2]). It is known for 
positive couplings J that the formation of ferromagnetic 
order needs a distinct coupling strength JS. For low cou- 
plings JS ferromagnetism only occurs at low band fillings 
n. The higher the coupling the higher the band filling 
can be before the ferromagnetic order breaks down. A 
maximum point is reached at half-filling n = 1 where fer- 
romagnetism never existsi^ii^i^i"— . Thus before we can 
compare different configurations of coupled films we have 
to find out if the intrinsic ferromagnetic order is stable. 
This is done in this work by a comparison of the ener- 
gies of the ferromagnetic {M = S) and the paramagnetic 
{M = 0) solutions. Normally at n « 1 antiferromagnetic 
phases occur which we neglect in this work. 
In uncoupled FM films neither parallel nor antiparallel 
alignment of magnetizations is favored, since the total 
QDOSs (i.e. spin-up plus spin-down) and with it the in- 
ternal energies are the same. Things change in coupled 
systems. The layers interact with each other via a hop- 
ping ej^(k||) in (|12p perpendicular to the film plane. This 
leads to a hybridization between the single layers' den- 
sities of states, resulting in a transfer and a repulsion of 
states to other energies (cf. Fig. [3]). The strength of the 
hybridization, i.e. the number of transferred states and 
the amount of repulsion, depends on the original (ener- 
getic) distance of the states of the uncoupled films - the 
closer the states the stronger the hybridization. 
Especially parallel and antiparallel configurations differ 
in coupled layer systems, because the hybridization acts 
only between the QDOSs of the same spin. For the par- 
allel alignment the QDOSs per spin direction of the two 
films are energetically close to each other, while they are 
far away for the antiparallel configuration. That is why 
the QDOS of the system with parallel alignment is con- 
siderably changed due to the hybridization. This can be 
directly seen for the two layer case in Fig. |31 

Changes in the QDOS are also changes in the inter- 
nal energy. So one can expect a strong influence of the 
hybridization on the phase diagram. We consider three 
types of phases: 

• both FM layers are ferromagnetically saturated 
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Figure 3: Quasi particle density of states of the top (FM), middle (NM) and bottom (FM) layer for parallel {solid line) and 
antiparallel configuration (dashed line). Differences to uncoupled layers occur due to a hybridization between the layers (cf. 
the magnified QDOS of the NM layer for Nl =3). At large couplings (shown here) the FM layers are energetically far away 
from the NM layers and the hybridization is very small. A special case is Nl ~ 2 where no NM layer exists and the FM layers 
hybridize directly. Vertical lines mark the chemical potential. Parameters: S = 7.5, JS — 2eV, n — 0.5, W^^ = leV, Tq^ — 



(M — S) and both magnetizations are aligned par- 
allel to each other 

as before but with antiparallel alignment 

the FM layers are not magnetically ordered (M 
0, paramagnetic 



Indeed all three phases occur in the phase diagram (cf. 
Fig. S]). Let us first discuss the Nl = 2 system, where 
the total system consists of two FM layers only. As in 
a bulk syste m^-'^i^^" — parallel alignment (corresponding 
to the ferromagnetic phase in the bulk) exists always at 
low doping n and increases its density range for strong 
couplings J. At n < 1 the antiparallel alignment is pre- 
ferred which can be compared to layered antiferromag- 
netic phases in the bulk^. Thus the parameter depen- 
dencies on J and n stay qualitatively the same as in the 
bulk system. As said before not all phases which have 



been found in bulk systems are considered in this work 
for simplicity. 

When we go to > 2 the FM layers are separated by 
Nl — 2 NM layers. Two new effects occur in this case. 
First it can be seen in Fig. |3] that for strong couplings 
parts of the QDOS of the FM layers are far below the 
QDOS of the NM layers. This means that the FM lat- 
tice sites are occupied by electrons at first which leads to 
a charge transfer from the NM to the FM layers. Thus 
the total electron density n normally differs from the den- 
sity within the FM layers. The total density in the 
system is 



(54) 



where A^°' is the absolute number of electrons and N is 
the total number of lattice sites. At strong couplings only 
the FM sites are occupied leading to a density in the FM 
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Figure 4: Phase diagram of the orientation of the magnetizations of the ferromagnetic layers with Nl ~ 2, 3, 4, 5. The 
magnetizations are either parallel (white) or antiparallel (striped). Additionally, areas where the layers themselves are not 
magnetically ordered are shown (gray). Vertical lines mark the critical band occupation where the FM layers are completely 
filled at large couplings (no = 2/Nl). Parameters; S = 7.5, W'^° = leV, T^^ = 



layers 



FM 



(55) 



with the number of lattice sites of the FM layers N^^. 
In our sandwiched layer systems only two FM layers 
are present in an Nl layer system which means that 
N™/N = 2/Nl. This means that the FM layers are 
half-filled 



2/Nl 



nc 



(56) 



Indeed this band occupation nc marks the transition 
to the paramagnetic region at large couplings. A sim- 
ilar effect has been found for the diluted Kondo-latticc 
model^. 

The second effect is due to the hybridization. When more 
layers are present in a system the hybridization leads to 
more complex shapes of the QDOS (cf. for example the 
middle layer in Fig. This manifests itself by a more 



complex structure of the phase diagram at low couplings 
or n < nc. Areas of parallel and antiparallel configu- 
rations of systems with larger Nl alternate more often 
than in systems with small Nl. Especially the increased 
number of peaks in the QDOS, compared to independent 
layers, influences the internal energy and with it the pre- 
ferred configuration. 

But there is another question. Even though we get 
a larger complexity of the QDOS with increasing Nl 
the absolute differences between the QDOSs for the two 
alignment types get smaller (cf. Fig. |3l for Nl = 5 the 
total QDOSs (spin-up plus spin-down) of the two align- 
ments are almost equal). From the phase diagram we 
do not know how large the energetic differences between 
these alignments are. Thus we define the interlayer ex- 
change 



AC/ = U^^ - U^^ 



(57) 
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Figure 5: (color online) Interlayer exchange A(7 — U^^ — U^^ 
for Nl ~ 2 (top), 3 (middle) and 4 (bottom) layers in depen- 
dence of carrier density n and coupling JS. Negative values 
(red) mean ferromagnetic and positive (blue) antiferromag- 
netic exchange. Parameters: S = 7.5, W^^ = leV, M = S, 



charge transfer from NM to FM layers. This results in 
a favored occupation of the FM layer and leads to half- 
filling at nc = 2/Nl. This changes when the NM QDOS 
lies a priori at low energies due to < 0. The influ- 

ence of a different position on the QDOS and the charge 
transfer can be seen in Fig. [51 When the NM QDOS 
comes in the vicinity of the FM QDOS it also gets occu- 
pied, leading to a lower density n^^. Especially a half- 
filling can be removed. Thus even for n > 2/Nl and large 
couplings ferromagnetically ordered phases can appear 
(cf. Fig. El). We use the special choice T^^ ^ -\JS 
which means that the NM QDOS is always at the lower 
subband of the FM QDOS. Thus the charge transfer does 
not change at strong couplings any more. Then the sys- 
tem gets in the usual double exchange limit, which means 
that the phase boundaries do not vary with J . However 
the number of phase areas of parallel and antiparallel 
alignment still increases with the number of layers ■ 
Not only the charge transfer influences the phase diagram 
but also the hybridization strength. With increasing cou- 
pling and fixed a position of the NM states near i? w the 
NM/FM layers' QDOSs moved apart from each other and 
the hybridization was reduced. This especially infiuenced 
the magnitude of the interlayer exchange AU . When we 
again use T^^ = -jJS the NM and FM QDOSs al- 
ways stay in the same distance and the hybridization 
strength does not change. So the interlayer exchange is 
not damped with increasing coupling any more (cf. Fig. 
[7]). Instead of that it increases its magnitude at low cou- 
pling and saturates at strong couplings. The oscillatory 
behavior with respect to the electron density is still pre- 
served. 



B. Curie temperatures 



as the difference between the two energies of both con- 
figurations. In this definition a negative AJ7 means a 
ferromagnetic interlayer exchange. We indeed see the ex- 
pected dependencies in Fig. [S] Even though an intrinsic 
magnetic order of a single FM layer is enhanced^^ by a 
large J, the interlayer exchange is heavily damped. This 
is due to the reduction of hybridization effects which are 
responsible for the intralayer exchange. For high Nl the 
damping is larger because an electron has to cross more 
NM layer to get from one FM layer to the other. The 
interlayer exchange has its maximum (absolute) value 
at intermediate couplings where hybridization effects are 
still relatively strong and the intrinsic ferromagnetism is 
already comparatively stable. 

Of course it is not required that the FM and NM layers 
have the same band's center of gravity Tq. In fact these 
two layer types consist in reality of different materials 
which have another band structure. Thus it is interest- 
ing to have different relative positions of the layers. We 
will focus on the choice = and an arbitrary JJNM 
in this work. 

One of the major effects at strong couplings J is the 



It is somehow contradictory to calculate transition 
temperatures in dimension reduced systems, i.e. in 
Z? < 3. The well-known Mermin- Wagner theorent^S tells 
us that there can be no ordered magnetism in films in 
the Heisenberg model at finite temperature. This also 
holds for some other many-body models and especially 
for the Kondo-lattice model^i^. But the statements 
are based on the existence of highly symmetric systems, 
which are not available in reality. Several symmetry- 
breaking effects (e.g. spin-orbit coupling—) can occur in 
the experiments or theory which results in finite tran- 
sition temperatures^"—. In fact we assume a break of 
symmetry in the derivation of the entropy Sq (M) where 
we assume that states with different magnetizations are 
non-degenerated. This can be done for example by a (vir- 
tual) magnetic field (cf. Sec. IIIB 2p . Additionally there 
can be NM layers in the total system which break the 
symmetry of the pure KLM, too. Thus the existence of 
Tc > in these films is indeed possible. 
It has been seen in the previous section that the param- 
eters J and n have a large influence on the stability of a 
parallel or antiparallel alignment of the magnetizations in 
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t7 = -JS/4 



t7 = -JS/2 




-10 1 2-2-10 1 2-2-10 1 2 
E[eV] E[eV] E[eV] 

Figure 6: Quasi particle density of states of the middle (NM) and bottom (FM) layer for parallel {solid line) and antiparallel 
configuration (dashed line) and different positions Tq^ of the NM layer. Vertical lines mark the chemical potential and the 
according carrier densities of the FM layers are shown. Parameters: S = 7.5, JS = 2eV, n = 0.8, W^^ = leV, Nl = 3 



the FM layers to each other. We now want to investigate 
the behavior at finite temperatures. To do this we have 
to use the free energy as the thermodynamical potential. 
In our work it was given as (|33]) 

Fm{T) ^Um{0) - T(5m(0) + /m(T)) (58) 

which can be mapped on a magnetization dependent 
function 



Ft{M) =Uo{M) ~ T{So{M) + It{M)) 
With the condition 

dFT{M) 



dM 







we get 



T(M) 



dM'Uo{M') 



dM'So{M')+dM'lT{M') 



(59) 



(60) 



(61) 



M 



This is the temperature when the system has the magne- 
tization M. The temperature dependent term dMlxiM) 
prevents a direct solution, but can be calculated in prin- 
ciple. It turns out in the Kondo-lattice model that it 
always holds dMSo{M) » dAjIriM) and we can neglect 



dAiIriM) in a very good approximation. This simplifies 
to 



T{M) = 



dM'UoiM') 



dM-So{M') 



(62) 



M 



where the right-hand side has no temperature depen- 
dence any more. Thus we get a direct equation for 
the temperature and with it for the magnetization curve 
M{T) = T~^{M). Especially we get the Curie tempera- 
ture by choosing 



dApUo(M') 



dM'SoiM') 



(63) 



M=0+ 



Let us discuss principal aspects of the formula. At high 
spins the entropy of the localized moments is much larger 
than that of the electrons. So we can set So{AI) = 
Sl°%M) + 5'§'(M) w S'i°"(M). The entropy of the local 
moments depends only on the magnitude of the magneti- 
zation \M\ and the spin quantum number S. This means 
that there do not occur any differences between paral- 
lel and antiparallel alignment in the entropy. Secondly 
the differences of the internal energy have to be taken at 
M ?a 0. At M = there is no difference between the 
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Figure 7: (color online) left: Phase diagram as in Fig. |4]but with Tq^ = —^JS for Nl ~ 3 (top) and 4 (bottom) layers, right: 
Interlayer exchange coupling as in Fig. [5]but with Tq^ = —\JS for Nl = 3 (top) and 4 [bottom) layers. Negative values {red) 
mean ferromagnetic and positive (blue) antiferromagnetic exchange. 



configurations and with it between the energies. It holds 

AmC/o"(M)|m=o+ =C/o"(AM > 0) - U^iM = 0) (64) 
=C/o"(AM > 0) - C/o(M = 0) (65) 

with a =tt7 ti- Thus we can calculate the difference of 
the Tcs to 

ATc = - (66) 

^ [/o^^(AA/>0)-t/j^(AM >0) 

AM'5o(Af')|M=o+ ■ ^ ^ 

The numerator of this fraction is the interlayer exchange 
AU defined in dST]) and ATc ^ ^U. The only difference 
is that we calculated AC/ in Sec. IIII Al at AI ~ S while 
in this case we use a small but finite AI — AM > 0. The 
entropy in ([S7)) plays more or less the role of a factor of 
proportionality. 

It is well known from RKKY calculation o^'^^ and 
experiments^ that the interlayer exchange coupling usu- 
ally has an oscillatory and damped behavior according to 
the spacer thickness. The RKKY calculations are based 
on the Heisenberg model which we do not use in our work. 
However also in our treatment the oscillations of the in- 
terlayer exchange are clearly visible (Fig. [S]). In fact the 
damping goes approximately with the square of the in- 
verse distance as in the RKKY treatment". Only 



the direct interlayer exchange between the two FM films 
(Nl — 2) is significantly larger than the others. To our 
opinion the RKKY behavior at Nl > 2 is recovered be- 
cause the RKKY is valid for small couplings J and in the 
PM layers J equals zero. There exists only an indirect 
interaction due to the hybridization with the FM layers 
which is small, too. Thus the conditions for an RKKY 
behavior are fulfilled. In contrast to that a direct cou- 
pling of two FM layers at large J is beyond the validity 
of RKKY calculations. Indeed calculations with modi- 
fied RKKY integrals show that the interactions range is 
significantly shorter at strong couplings^^. 
The same oscillations as for AU occur for the (relative) 

ATc = (T^^ -TI^)/{T^^ + T^^) (Fig. ^ which is closely 
related to the interlayer exchange as shown in (p7)) . The 
connection of the interlayer exchange to an oscillatory 
behavior of Tc has been found in other work^^ where the 
origin of this behavior is traced back to quantum-well 
states. As for the case of AU at vanishing temperature 
the differences between T^tt and T^" are heavily damped 
if the NM layers are at T^^ = 0. 

When the NM QDOS is near the FM QDOS due to the 
choice Tq^'^ = —\JS the amplitudes of the oscillations 
remain relatively strong even for a large spacer thick- 
ness, i.e. a large Nl in our work. The reason is the same 
as for T = 0. Since the QDOSs are near to each other 
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Figure 8: (color online) left: Relative difference of Curie temperatures ATc = {Tq — T^y ) / {T^ +Tq) vs. number of layers Nl 
for T^^ = and T^^ = -i JS. Oscillations are much stronger when the NM layer is near the FM layer (i.e. T^^ = -i JS). 
Additionally the interlayer exchange dJJd? with the distance d — Nl — 1 for n = 0.9 is plotted (dashed line). For Nl > 2 it 
shows the typical RKKY-damping (AC/ ~ d~^). top right: Reduced (without It{M)) free energy Ft{M) = Uo{M) - TSo{M) 
for parallel {solid Ime) and antiparallel (dashed) alignment at n = 0.2. The antiparallel alignment is always lower or equal 
(M = 0) in energy. The optimal magnetization is at the minimum of Ft{M) which is for the shown curves at M°p' = or 
M°P' _ ^ except for antiparallel alignment at T = 20K (M"^' = 5.1, see arrow), bottom right: Resulting free energy at the 
optimal magnetizations. The alignment with the higher Curie temperature (here T^'^ ~ 40K > T^t^ « 15K) is always lower 
(T < T^,^) or equal (T > T^^) in energy. Parameters: S = 7.5, W^^ = leV, JS = leV 



the hybridization is strong and there are large differences 
between both alignment types. This leads to a large in- 
terlayer exchange and with it to a large difference in the 
Curie temperatures. 

It has to be said that often both configurations of the 
magnetizations of the films have a positive transition 
temperature. Thus both can be the energetically opti- 
mal phase at finite temperatures and temperature driven 
transitions between both alignments could be in principle 
possible. But usually the phase with the higher transition 
temperature exists in the whole temperature range 
[0,r^]. This can be shown by simple considerations. It 
has been stated before that there is no difference between 
the configurations at M = 0. Thus also the energies of 
the both types of alignment are equal to each other 

U^'^iM = 0) = U^^{M = 0) . (68) 

Secondly Uo{M) and also dMUQ{M) are monotonic func- 
tions of M (cf. Fig [S]) which means that the phase with 
the lowest internal energy at ferromagnetic saturation 



M — S also has the lowest energy for all M. Due to 
the fact that the entropy So{M) and the integral It{M) 
in ([59)1 hardly depend on the type of alignment the free 
energy is lower for all M > and T and the transition 
temperature is higher, too. Thus a phase with higher Tc 
has a lower free energy for all A/ > and is preferred for 
all T (cf. Fig. HI). 

Since the differences of the Curie temperatures are 
strongly connected to the interlayer exchange, all pa- 
rameter dependencies are similar to the dependencies of 
AC/. The most influential parameters in the Kondo lat- 
tice model are the coupling J and the band occupation 
n. 

Generally in bulk systems Tc increases with J for low 
and intermediate couplings and saturates at high cou- 
plings ("double exchange limit"). This is also true in 
the film systems treated in this work (cf. Fig. ^ 
whereas the saturation effect is more visible for the case 
of TqNM = -\JS. At T^^ = the carrier density n™ in 
the FM layers gets larger with increasing J due to charge 
transfer. For those higher densities the saturation tends 



13 




Figure 9: (color online) Curie temperatures in parallel {solid lines) and antiparallel {dashed) configuration for Tq^^ = {black) 
and Tcf^ = ~\JS {red, symbols), left: Tc vs. JS at n = 0.3. For T^^ = -\JS differences between the Tcs stay large even 
for strong couplings J. right: Tc vs. n at JS = leV. For — —^JS oscillations are stronger and a finite Tc exists for large 
n due to suppressed charge transfer. Parameters: S = 7.5, W^^ = leV 



to start at larger couplingsiii^. Again the difference be- 
tween the Curie temperatures of both configurations is 
interesting for coupled film systems. We see the typical 
behavior. For JJ^m _ g ^j^g difference decreases even 
though the absolute values of Tc increase. This is due to 
the reduction of the hybridization strength because of the 
larger distance between the FM and NM layer's QDOS. 
For Tq^ — -~\JS the hybridization remains approxi- 
mately the same and therefore the (absolute) difference 
between both Tps increases with raising transition tem- 
peratures. 

The role of the charge transfer is more obvious when we 
look on the influence of n on Tc- It can be seen in Fig. [9] 
that the range of finite Tc according to n is reduced with 
Nl when we choose T^^ — 0. It has the same reason as 
in the case of T = where only paramagnetism occurred 
at large couplings and n > nc = due to the half 

filling of the FM layers (cf. Fig. HJ. Since charge transfer 
to the FM layers is reduced for Tq^ = —^JS or rather 
there is even transfer to the NM films, the density in the 
FM layers stays at a moderate level < 1 even at 

large n. Thus also Tc remains positive over the whole 
range of n £ [0,1] for intermediate and strong couplings 
J. Again the maximum differences between and T^^ 
remain large at a high layer number Nl- As in the case 
of r = for the interlayer exchange AU the Tcs are 



strongly oscillating with changing n. 



IV. CONCLUSIONS 

In this work we investigated thin film systems with re- 
spect to the alignment of the magnetizations of ferromag- 
netic layers. We concentrated on "sandwiched" system 
where two ferromagnetic (FM) layers are separated by 
an arbitrary number of non-magnetic (NM) layers. The 
magnetizations of the FM layers could either be parallel 
or antiparallel to each other. One main focus was on the 
influence of carriers on the preference of parallel or an- 
tiparallel alignment of the magnetizations. Thus we used 
a Kondo-lattice model to describe interactions of elec- 
trons with local moments within the FM layers. With an 
equation of motion method we could calculate the elec- 
trons' properties and with it the internal energy of the 
system at vanishing temperature T — 0. By a compar- 
ison of the internal energies of the parallel/antiparallel 
alignments we got magnetic phase diagrams in depen- 
dence of the carrier density n, the coupling J and the 
number of layers N^. With increasing number of layers 
there are more alternations of the stable alignment ac- 
cording to n. A very important quantity is the interlayer 
exchange which we defined as the difference of the inter- 
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nal energies AI7 = [7^^ — /7^~'-. Depending on the position 
of the NM layers the interlayer exchange can either 
be constant or heavily damped according to an increase 
of the coupling J. 

Additionally we presented a method to calculate a mag- 
netization dependent free energy Ft{M) without any use 
of a Heisenberg model. When we minimize Ft{M) ac- 
cording to M we find the optimal magnetization for any 
given temperature T. With it we can calculate Curie 
temperatures T^^, tJ}' of the respective alignments of the 
magnetizations. These, or better the differences between 
them, are strongly connected to the interlayer exchange 
AJ7. As for AC/ the differences between the Tqs show 
an oscillating behavior according to n. The number of 
oscillations depend on the number NM layers. 
The results in this work show, that the carrier density 
heavily influences the interlayer exchange in the Kondo 
lattice model, which is maybe not captured in the stan- 
dard RKKY Heisenberg model treatment. Ab-initio cal- 
culations of real materials are in principle exact at T = 
and contain more interactions, e.g. superexchange, which 
makes a direct connection to the pure KLM difficult. But 
also in ab-initio calculations a mapping onto model sys- 
tems, e.g. the Heisenberg model, is necessary to get re- 
sults for T > 0. This can change the general behavior 
of the system being the reason why we wanted to avoid 
such a mapping. 

It would be very interesting to extend the number of FM 
layers in the investigated systems or to add semi-infinite 
leads to the top and bottom layers. This would allow a 
comparison to bulk systems or potentially lead to new 
effects. Since these calculations would be numerically 
extensive they are left for future work. 

Appendix A: Differences to Mean-field results 

The simple mean-field approximation (MFA) of the 
KLM is often used in literature. In this approximation 
the spin-fiip terms of the Hamiltonian ([T]) are neglected 
resulting in a self-energy 

Mf^'^E) - -z^ijM . (Al) 

At T = 0, or more precisely at ferromagnetic satura- 
tion M = S, the spin-flip part plays a minor role. Thus 
the T — Q phase diagrams remain almost unchanged in 
the MFA except that the paramagnetic phase is missing. 
This has also be found for the pure and diluted KLM^^i^^. 
When the magnetization is away from saturation the 
spin- flip part gets important. This is especially true at 
Tc where the magnetization vanishes and it can be ex- 
pected that the MFA has several flaws. As an example 



we want to discuss the J-dependence of Tc (cf. Fig. [TOl) . 
Indeed the results are different from those of Fig. |9l For 
T^^ = we get the typical MFA dependence Tc 
which means we get no saturation of Tc at large J. Sec- 
ondly, and maybe more important in terms of interlayer 
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Figure 10: (color online) Curie temperatures in the mean- 
field approximation for parallel (solid lines) and antiparal- 
lel (dashed) configuration for T^^ = (black) and T^^'^ = 
— I JS* (red, symbols). The results differ drastically from those 
of Fig. H Parameters: S = 7.5, W^'^ = leV, n = 0.3, Nl ^ 3 



exchange coupling, the difference of the Curie tempera- 
tures of different alignments increases with increasing J. 
For T^^ = this behavior is reversed - ATc de- 

creases with increasing J. This are completely contrary 
results compared to the MCDA calculations. 
The reason for that can be easily understood. Due to the 
simple structure of the self-energy of the MFA the sub- 
bands of the QDOS are at -E w —Za^JM which means 
their positions depend on the magnetization. Especially 
for M = at Tc the sub-bands are at E; w 0. The differ- 
ence of the Tcs depends on the interlayer exchange which 
is large when the QDOSs of the FM and PM layers are 
close to each other. In the MFA this is always the case 
for M = and T^^ = 0. Thus ATc is large for every J 
and its absolute value increases with the absolute value 
of Tc. For = -i J5 the QDOS of the PM layers 

can be far away from the FM layers' QDOS for a large 
J and thus ATc decreases. Because in the MCDA the 
subbands of the FM layers stay for aU M at i? w ±^JS 
the situation is inverted - the FM/PM QDOSs are close 
to each other for T^'^ — ^JS and far away from each 
other for T^^ = 0. 

Thus, besides the drastic overestimation of Tc i the mean- 
field approximation can lead to completely wrong predic- 
tions in the KLMi^ and should be handled with care. 
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